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Abstract 

We examine capabilities of the modified perturbation theory (MPT) for de- 
scription of the processes with productions and decays of fundamental unstable 
1 particles. We calculate total cross-section for e + e~ — > 7, Z — > W + W~ — > 4/ 

in a model with the Dyson-resummed and with the MPT-expanded up to the 
NNLO Breit-Wigner factors, and compare the outcomes. At the ILC energies a 
Q_i' coincidence of the outcomes is detected with precision better than one per-mille. 
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An investigation of the W^-pair production plays an important role in testing the 



Standard Model and searching for physics beyond. At the International e + e~ Linear 
Collider (ILC) pQ the accuracy of the appropriate measurements is planned at the 
■ per-mille level. This implies that the theoretic calculations of the cross-section for W- 
. pair production must be carried out with the per-mille precision or better. To achieve 
this, the calculations in the next-to-next-to-leading order (NNLO) are required and, 
in addition, the resonant contributions of unstable particles must be treated with the 
appropriate accuracy without violation of the gauge cancellations. 

A simultaneous implementation of the above requirements is a very difficult task, 
unresolved so far. Really, the approach of the double pole approximation (DPA) 
successfully applied at LEP2, provides only the NLO precision and in the resonant 
region only [2]. So the DPA is unable to solve the problem [3]. The pinch-technique 
method and the background- field formalism, another for a long time developed ap- 
proaches [4} [5] , seemingly can provide the NNLO precision. But for the maintenance 
of the gauge cancellations, they require a lot of calculations of extra contributions 
that appear formally beyond the necessary precision, which is impractical [6]. So, at 
present the main hopes are pinned on the approach of "complex-mass scheme" (CMS) 
which avoids the mentioned difficulties [31 [7]. Nevertheless, the CMS due to intrinsic 
problem related to unitarity, provides the NLO precision only [7J. So, alternative 
approaches are required. 

A modified perturbation theory (MPT) [HI El [10] is a probable candidate to become 
such an approach. Its main feature is the direct expansion of the probability instead 
of amplitude in powers of the coupling constant with the aid of distribution-theory 
methods. The latter methods allow one to impart a well-defined meaning to the 
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resonant contributions in the expansion of probability. A condition of asymptoticity 
(and completeness) of the expansion must ensure the gauge cancellations. In the case 
of pair production of unstable particles, the most-elaborated description of the method 
was given in [TU] . In [TTj [T2"] the convergence properties of the MPT expansion were 
tested in a model related to the top-quark pair production. At the ILC energies, a 
good convergence of the MPT series was detected with the precision of the description 
within the NNLO up to a few per-mille. 

In this Letter we present results of similar investigation in a model for VT-pair 
production in e + e~ annihilation. We emphasize that a separate analysis in this case 
anyway must be made since this process is planned for the ILC, and since this pro- 
cess has peculiar features. Namely, in the case of VT-pair production there appear 
the i-channel contributions, which in the cross-section result in the contributions 
with a more complicated analytic structure (the additional logarithm, see e.g. [13J). 
Moreover, there appear the large gauge cancellations arising even in the Born approx- 
imation from contributions of longitudinally polarized P^-bosons. In particular, the 
mentioned cancellations reach one order of magnitude at 400 GeV and two orders at 

1 TeV [13] . So, one should check the behavior of the MPT in the presence of such 
large cancellations. 

In the next section, we detail the model in which we carry out calculations. The 
numerical outcomes are presented in Sec. [31 In Sec. HI we discuss the results. 

2 The total cross-section and a model for VF-pair 
production 

The total cross-section of H^-pair production in e + e~ annihilation has the form of a 
convolution of the hard-scattering cross-section with the flux function [T3| [14"]. 



(The angular distributions may be described in the MPT, as well, but we do not 
discuss this option here.) The flux function stands for contributions of nonregis- 
tered photons emitted in the initial state. Below we consider in the leading-log 
approximation (see details e.g. in [13]). The hard-scattering cross-section is given by 
an integral over the virtualities of unstable particles, 



a(s) = / / dsids 2 0{y/s-y/si-y/s2)\/\(s,Si,S2) $(s;si,s 2 ) (1+8 C ) p(sx)p(s 2 ) ■ (2) 



Here s lmin and s 2min are minima of the virtualities, y/su^ + y/s 2m i n = ^/s min . The 
A = [s — (y/si +a/s2) 2 ][s — (y/si — v^2 ) 2 ] is the kinematic function. The ^-function 
and the square root of A constitute a kinematic factor. The p(sj) are Breit-Wigner 
(BW) factors of W^-bosons. The <3> and the factor (1 + S C ) stand for the amplitude 
squared and soft photon contributions. Since $ corresponds to one-particle irreducible 
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contributions, it has no singularities on the mass-shell of H^-bosons. On the contrary, 
if we naively expand the BW factors in powers of the coupling constant, nonintegrable 
singularities will be generated. This fact constitutes a well-known problem of the 
determination of the cross-section in the case of productions and decays of unstable 
particles. In the framework of systematic expansion of amplitude in powers of the 
coupling constant this problem, seemingly, is unsolvableQ However, at the level of 
the cross-section, a solution to the problem does exist. 

The basic idea is to consider the singularities in the cross-section in the sense 
of distributions. In a mathematically correct way, the problem may be stated as a 
problem of asymptotic expansion of BW factors in powers of the coupling constant in 
the sense of distributions. In the case of a separately taken BW factor with smooth 
weight, a solution to the problem is well-known [8]. Namely, the expansion is begin- 
ning with the 5-function, which corresponds to the narrow-width approximation. The 
naive Taylor-expansion terms are supplied with the principal-value prescription for 
the poles. Strongly nontrivial contribution appear as a sum of the delta-function and 
its derivatives with coefficients c n , which are polynomials in the coupling constant a. 
Within the NNLO, the expansion looks as follows: 
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2aReSi(s) 



(s-M 2 ) 2 (s-M 2 f 



(3) 



+ y c n M- 5 {n) is-M 2 ) + 0(« 3 ) . 

Here M is the renormalized mass of the unstable particle, To is its Born width, 
E = a Ei + a 2 E 2 + a 3 E 3 + ■ • • is the self-energy in the scalarized propagator. Co- 
efficients c n within the NNLO are determined by three-loop self-energy contributions 
and by their derivatives, all taken on-shell. The structure of the contributions is such 
that in the OMS-type schemes of the UV renormalization the real self-energy contri- 
butions appear without the derivatives or with the first-order derivative only. Such 
contributions are determined by the renormalization conditions. In the unstable- 
particles case a convenient scheme of the UV renormalization is the OMS or "pole" 
scheme [151 HE], defined by equating the renormalized mass of unstable particle to real 



component of the pole of its propagator. This implies that the OMS mass coincides 
with the observable mass, which is independent of the UV renormalization scheme 
(and which is gauge- invariant). The second renormalization condition is defined by 
equating imaginary component of the on-shell self-energy to imaginary component of 
the pole. Coefficients c n in this scheme are determined as follows [TU] : 



12 , 2 
Co = — a — + a 

h 



r IJ_h 

n h 



in 



f\2 



ci = 0, 



c 2 



^2 t-2 



(4) 



1 Let us remind that the CMS cannot be considered as a systematic approach for the calculation 
of amplitude beyond the NLO [7] . 



3 



where I k = ImE fe (M 2 ), and I{ = lmE^(M 2 ). 

The next ingredient of the MPT is the analytic regularization of the kinematic 
factor via the substitution [A(s, Si, s 2 )] 1 ^ 2 — > [A(s, si, s 2 )] 1/ [ID]. The regularization 
is necessary to impart enough smoothness to the weight at the BW factors in for- 
mula 02]). Under the regularization, by representing the remaining contribution to 
the weight in the form of Taylor expansion with a remainder, we can do analytic 
calculation of singular integrals irrespective of details of the definition of the weight. 
Further, by putting v = 1/2, we obtain finite expressions with the expansion being 
asymptotic [10J. Thereby the problem is reduced to numerical calculations only. 

Now we turn to the definition of the model for carrying out the calculations. 
Acrually, we define the model in many ways by analogy with [12]. First of all, we note 
that the MPT-expansion of the BW factor based on the full propagator of PU-bosons 
coincides within the NNLO with that based on the following modelling propagator: 

A _1 (s) = s - M 2 + aReSi(s) + ialmEi(s) 

+ a 2 [R 2 + H 2 + (s-M 2 )(R' 2 + H' 2 )] +a 3 (R 3 + H 3 ). (5) 

Here R k = ReEfc(M 2 ) and R' k = ReE-[(M 2 ). When considering in the sense of con- 
ventional functions, a discrepancy between propagator (jSJ) and the full propagator 
turns out to be beyond the NNLO, as well. Below we consider propagator (jSJ) as the 
full or "exact" propagator in the model. The ReE 1 (s) and ImEi(s), we define by 
means of direct calculations in the SM. However, for sake of simplicity we restrict 
ourselves by consideration of only quark and lepton contributions. In this way we 
avoid the IR divergences generally arising at calculating ReEi(s). For the definition 
of the on-shell self-energy contributions in <^j, we take advantage of the UV renor- 
malization conditions (we use the OMS scheme) and the conditions of unitarity. In 
this way we obtain (see details in [12] and [15]) 

R 2 = -hl[ , R' 2 = -hI'{/2 , R 3 = -I 2 I[ - I X I' 2 + IfR>;/2 , (6) 
aI 1 = MT , a 2 I 2 = MT 1 , a 3 I 3 = MT 2 + T 3 /(8M) , (7) 

where r , r 1; T 2 are the Born-, one-loop, and two-loop contributions to the width. 
The I 2 , we determine with the aid of an approximate relation a 2 I 2 = (Tx/Tq) al[. 

Further we proceed to the definition of the amplitude squared $ in formula (J2]). 
For the purpose of testing only of the convergence properties of the MPT expansion, 
we determine $ in the Born approximation. In the main we follow [13] . except Z- 
boson propagator. Recall that in the Born approximation the Z-boson propagator 
arises immediately owing to e + e~ annihilation, and that it makes one-particle re- 
ducible contribution. Therefore this contribution becomes resonant if s' in formula 
([I]) falls on s' = M§ . In the latter case the mentioned propagator of Z-boson in $ 
is to be MPT expanded]! Unfortunately, the loop corrections to Z-boson propagator 

2 Let us remark that the corrections to $ outside the above mentioned Z-boson propagator, 
constitute one-particle irreducible contributions. Therefore they are non-resonant and should be 
considered by means of conventional perturbation theory. 
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without simultaneous consideration of the radiative corrections to vertices, will lead 
to violation of unitarity at high energies. For this reason, we restrict the energy range 
in which we consider corrections to Z-propagator. Namely, we determine Z-boson 
propagator by formulae ©-dZj) at s < sz, and as a free propagator at s > sz, where 
sz is a parameter such that sz > Mf . By means of this trick we, in a proper way, 
take into consideration the resonant contributions of Z-boson propagator, and simul- 
taneously preserve the gauge cancellations and the unitarity at high energies. (Let 
us remember that this is a model consideration. In the approach with the taking 
into consideration of all radiative corrections, both in the vertices and in the Z-boson 
propagator, the gauge cancellations and, consequently, the unitarity should automat- 
ically be maintained. This should be so, since we construct a complete asymptotic 
expansion in powers of the coupling constant.) 

For completing the definition of the model, we must determine also the soft- 
photon factor (1 + S c ) in formula (121) . Let us remember that we have omitted photon 
contributions in the self-energies of H^-bosons. So we have to ignore those soft- 
photon contributions whose IR-divergent contributions are to be cancelled in the 
cross-section. For this reason, we consider only Coulomb soft-photon contributions in 
5 C . Concretely, we take advantage of the well-known formula for the Coulomb factor in 
the one-photon approximation, which includes appropriate off-shell and finite-width 
effects, see details in [TO] . 

Now the model is completely determined, and the cross-section in the model 
may be straightforwardly calculated. We call the outcome the "exact result in the 
model". After the MPT-expansion, the cross-section is represented in the form 
er(s) = o- LO {s) + aa\(s) + a 2 cr 2 (s), where a^o is the cross-section in the LO ap- 
proximation, a a i and a 2 o"2 are the NLO and NNLO corrections, respectively. So, 
the onlo = °~lo + o. o"i and onnlo = c"lo + a<7i + o^oi determine the NLO and the 
NNLO approximations. 

3 Numerical results 

For numerical calculations, we use parameters: M w = 80.40 GeV, M z = 91.19 GeV, 
rrit = 175 GeV, and the other quarks and leptons are massless. To parametrize the 
amplitude squared we use effective coupling determined as 

a G ^ = V2G^s 2 w /7r, (8) 

where s 2 w = 1-M^/Mf, is the Fermi constant, = 1.16637 x 10" 5 GeV^ 2 . For 
the calculation of relative corrections, we use a = 1/137.036. The parametrization of 
this kind was used e.g. in [7]. Since we consider the widths as the relative corrections, 
the mentioned parametrization implies large electroweak (EW) corrections to the 
widths. To estimate specifically the 0(a) EW corrections, we use results of [14] and 
|18j . The 0(a s ) corrections, we obtain in the "naive" approach, extracting them from 
the corresponding factor (1 + a s /ir) with a s = 0.12. In view of the lack of complete 
results on the two-loop corrections to the widths, we estimate they as the differences 
between the current values T w = 2.082 GeV, T z = 2.534 GeV [19] and the sums of 
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the appropriate Born and one-loop results 



In this way we obtain 



17 = 1.977 GeV, 
rf = 0.102 GeV, 
rf = 0.006 GeV, 



rf = 2.362 GeV, 
rf = 0.169 GeV, 
rf = -0.036 GeV. 



(9) 



Another way to estimate the two-loop corrections is to consider they as equal to the 
two-loop QCD corrections extracted from the factor (1 + a s /^ + 1.409 a 2 /7r 2 ) [19]. 
In this way we obtain Tf = rf = 0.003 GeV. Below we use the latter option for 
checking stability of the outcomes, but as the main estimate we consider (Q. Recall 
that we use the widths for the determination of the "exact" propagator fl5]) and of 
the coefficients c n in the MPT expansion of the BW factors. Barameter sz in the 
definition of Z-boson propagator, we determine on the basis of the condition of the 
most smooth "sewing" solutions above and below sz- This gives J~s~z = 130-150 GeV, 
and the ultimate result is practically independent from the concrete choice. Finally, 
we put -y/sjmin = 30 GeV in formula (j2J). The calculations, we carry out on the basis 
of BORTRAN code with double precision written in accordance with results of [10J. 

The outcomes of the calculations are presented in Figs. [TH3] and in Table 1. In 
Fig. [U the thick curve shows the exact result for the total cross-section, the dotted and 
long-dashed curves show the cross-section in the LO approximation of the MFT and 
in the DBA, respectively. Recall that the LO in the MFT coincides with the narrow- 
width approximation. The DBA is constructed by substitution of $(s; M 2 , M 2 )(l + 
5 c (s;M 2 ,M 2 )) for $(s;si,s 2 )(l + S c (s; s u s 2 )) and s — M 2 + i T for A _1 (s) in the 
expression for the exact result. The curves for the NLO and NNLO are not presented 
in Fig. [U since they merge with the curve for the exact result at a given scale. 
Instead, they are shown in Figs. |2] and EJ where the percentages with respect to 
the exact result are presented. In Fig. |3] the results of Fig. [2] are repeated with 
greater scale on vertical axis. In Table 1 the outcomes are presented in numerical 
form at the characteristic energies at the ILC. In the last column the numbers in 
parenthesis represent uncertainties in the last digits arising due to computations. In 
other columns the uncertainties are omitted as they are beyond the precision of the 
presentation of data. The estimation of uncertainties is made in accordance with 
the algorithm described in [2U]. The results of checking stability of the outcomes 
with respect to the variation of the corrections to the widths (see explanation in the 
previous paragraph) are shown in Fig. HI 

The presented outcomes exhibit stable behavior of the NLO approximation and 
very stable behavior of the NNLO approximation in the energy region beginning with 
approximately 200 GeV. At the same time, the accuracy of the NNLO approximation 
in this energy region (more precisely, above 220 GeV) is less than one per-mille. These 
results remain valid at varying the input-data for the widths Tf and Tf , cf. Fig. [3] 
and Fig. HI The corresponding shift of the NNLO approximation with the changing 
of the widths constitutes a few of 10~ 4 in relative units. The concurrent shift of the 
NLO curve exactly follows the shift of the "exact" solution determined by propagator 
(J5]), and in relative units constitutes a few of 10 -3 . As for the DBA, it has rather 
non-stable behavior and at high energies it becomes openly bad. In particular, at 
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Figure 1: Total cross-section: the exact result in the model, the results in the LO in the 
MPT (the narrow- width approximation), and in the DPA. 
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Figure 2: Percentages with respect to "exact" cross-section calculated in the LO, NLO, 
and in the NNLO approximations in the MPT, and in the DPA. 
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Figure 3: Percentages to the exact cross-section, calculated in the NLO and NNLO approx- 
imations in the MPT. The horizontal line at 100.0% means the exact result. The notation 
for the curves is the same that in Fig. [2l 
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0.2 


15.258 


17.839 


15.175 


15.235(2) 




100% 


116.92% 


99.46% 


99.85(1)% 


0.5 


6.9355 


7.5657 


6.9294 


6.9342(7) 




100% 


109.09% 


99.91% 


99.98(1)% 


1 


2.8286 


2.9733 


2.8263 


2.8285(3) 




100% 


105.12% 


99.92% 


100.00(1)% 


3 


0.61023 


0.55733 


0.60625 


0.61026(6) 




100% 


91.33% 


99.35% 


100.00(1)% 


5 


0.391580 


0.24353 


0.31047 


0.31566(3) 




100% 


77.11% 


98.31% 


99.95(1)% 



Table 1: The results of the calculation of the total cross-section in pb and in % with respect 
to the exact result in the model. 
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Figure 4: The same as in Fig. [3] but with the other input for widths (see explanation in 
the text). 



the energies above 1.5 TeV the DPA becomes much worse even of the narrow width 
approximation (the LO approximation in the MPT). 

In the end of this section it is worth noticing that at the approaching to the thresh- 
old of PF-pair production, the accuracy and stability of MPT rapidly become worse. 
Actually this behavior was predicted in [10] and was observed at previous calculations 
[TT1 [12] . To prevent this difficulty another mode of the MPT near threshold should 
be applied, which includes Taylor expansion of a(s) not only in powers of a, but also 
in powers of the distance between s and the threshold JTU]. A study of the behavior 
of the MPT in this mode is a task for future investigations. 



4 Discussion and conclusions 

We have tested applicability of the MPT for description of JF-pair production in e + e~ 
annihilation. We have found that MPT provides stability of outcomes at the energies 
beginning with approximately 200 GeV. By means of model calculations we have 
shown that at the mentioned energies the NNLO approximation in the MPT provides 
better than one per-mille precision of the description of the total cross-section. 

The obtained precision exceeds that observed in the case of the top-quark pair 
production [12]. The difference is caused by the the lower values of the corrections 
to the width of VF-boson as compared to the corrections to the width of the top- 
quark (and to a lesser degree by the opposite sign of the corrections). Specifically, 
if we use the corrections to the width of the top-quark in place of those for the W- 
boson (in relative units) in the present calculations, the outcomes for the cross-section 
in relative units are obtained almost identical to those in the case thoroughly with 
the top-quark pair production. On the one hand, this confirms the applicability of 
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the MPT in the presence of the large gauge cancellations that occur in the case of 
the VT-pair production. On the other hand, this confirms the previous observation 
that the results of the MPT calculations in relative units depend rather weakly on 
the particular form of the test function (the amplitude squared without the resonant 
factors), but depend mainly on the structure of the resonant factors [12]. In total, the 
above observations mean that the modification of the test function due to turning on 
the loop corrections, in relative units should not lead to any significant modifications 
in the cross-section. So, our main results should remain in force in the presence of 
the loop corrections, in particular, the result about the better than one per-mille 
precision of the description of the cross-section. 

In summary, we conclude that the MPT is really a good candidate for the de- 
scription of VT-pair production at the ILC For a definitive conclusion, of course, a 
more comprehensive analysis must be made, in particular the analysis of the angu- 
lar distributions. A study of complete realistic processes, such as CC10, CC11 etc., 
in the framework of an appropriate model must be made, as well. The mentioned 
investigations will be carried out elsewhere. 
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